The Cystic Fibrosis Upper and Lower Airway Metagenome.

The microbial metagenome in cystic fibrosis (CF) airways was investigated by whole-genome shotgun sequencing of total DNA isolated from nasal lavage samples, oropharyngeal swabs, and induced sputum samples collected from 65 individuals with CF aged 7 to 50 years. Each patient harbored a personalized microbial metagenome unique in microbial load and composition, the exception being monocultures of the most common CF pathogens Staphylococcus aureus and Pseudomonas aeruginosa from patients with advanced lung disease. The sampling of the upper airways by nasal lavage uncovered the fungus Malassezia restricta and the bacterium Staphylococcus epidermidis as prominent species. Healthy and CF donors harbored qualitatively and quantitatively different spectra of commensal bacteria in their sputa, even in the absence of any typical CF pathogen. If P. aeruginosa, S. aureus, or Stenotrophomonas maltophilia belonged to the trio of the most abundant species in the CF sputum metagenome, common inhabitants of the respiratory tract of healthy subjects, i.e., Eubacterium sulci, Fusobacterium periodonticum, and Neisseria subflava, were present only in low numbers or not detectable. Random forest analysis identified the numerical ecological parameters of the bacterial community, such as Shannon and Simpson diversity, as the key parameters that globally distinguish sputum samples from CF and healthy donors. IMPORTANCE Cystic fibrosis (CF) is the most common life-limiting monogenetic disease in European populations and is caused by mutations in the CFTR gene. Chronic airway infections with opportunistic pathogens are the major morbidity that determines prognosis and quality of life in most people with CF. We examined the composition of the microbial communities of the oral cavity and upper and lower airways in CF patients across all age groups. From early on, the spectrum of commensals is different in health and CF. Later on, when the common CF pathogens take up residence in the lungs, we observed differential modes of depletion of the commensal microbiota in the presence of S. aureus, P. aeruginosa, S. maltophilia, or combinations thereof. It remains to be seen whether the implementation of lifelong CFTR (cystic fibrosis transmembrane conductance regulator) modulation will change the temporal evolution of the CF airway metagenome.

Competitive fitness of Pseudomonas aeruginosa isolates in human and murine precision-cut lung slices.

Chronic respiratory infections with the gram-negative bacterium Pseudomonas aeruginosa are an important co-morbidity for the quality of life and prognosis of people with cystic fibrosis (CF). Such long-term colonization, sometimes lasting up to several decades, represents a unique opportunity to investigate pathogen adaptation processes to the host. Our studies aimed to resolve if and to what extent the bacterial adaptation to the CF airways influences the fitness of the pathogen to grow and to persist in the lungs. Marker-free competitive fitness experiments of serial P. aeruginosa isolates differentiated by strain-specific SNPs, were performed with murine and human precision cut lung slices (PCLS). Serial P. aeruginosa isolates were selected from six mild and six severe CF patient courses, respectively. MPCLS or hPCLS were inoculated with a mixture of equal numbers of the serial isolates of one course. The temporal change of the composition of the bacterial community during competitive growth was quantified by multi-marker amplicon sequencing. Both ex vivo models displayed a strong separation of fitness traits between mild and severe courses. Whereas the earlier isolates dominated the competition in the severe courses, intermediate and late isolates commonly won the competition in the mild courses. The status of the CF lung disease rather than the bacterial genotype drives the adaptation of P. aeruginosa during chronic CF lung infection. This implies that the disease status of the lung habitat governed the adaptation of P. aeruginosa more strongly than the underlying bacterial clone-type and its genetic repertoire.

iPSC culture expansion selects against putatively actionable mutations in the mitochondrial genome.

Therapeutic application of induced pluripotent stem cell (iPSC) derivatives requires comprehensive assessment of the integrity of their nuclear and mitochondrial DNA (mtDNA) to exclude oncogenic potential and functional deficits. It is unknown, to which extent mtDNA variants originate from their parental cells or from de novo mutagenesis, and whether dynamics in heteroplasmy levels are caused by inter- and intracellular selection or genetic drift. Sequencing of mtDNA of 26 iPSC clones did not reveal evidence for de novo mutagenesis, or for any selection processes during reprogramming or differentiation. Culture expansion, however, selected against putatively actionable mtDNA mutations. Altogether, our findings point toward a scenario in which intracellular selection of mtDNA variants during culture expansion shapes the mutational landscape of the mitochondrial genome. Our results suggest that intercellular selection and genetic drift exert minor impact and that the bottleneck effect in context of the mtDNA genetic pool might have been overestimated.

Genome-Wide Methylation Mapping Using Nanopore Sequencing Technology Identifies Novel Tumor Suppressor Genes in Hepatocellular Carcinoma.

Downregulation of multiple tumor suppressor genes (TSGs) plays an important role in cancer formation. Recent evidence has accumulated that cancer progression involves genome-wide alteration of epigenetic modifications, which may cause downregulation of the tumor suppressor gene. Using hepatocellular carcinoma (HCC) as a system, we mapped 5-methylcytosine signal at a genome-wide scale using nanopore sequencing technology to identify novel TSGs. Integration of methylation data with gene transcription profile of regenerated liver and primary HCCs allowed us to identify 10 potential tumor suppressor gene candidates. Subsequent validation led us to focus on functionally characterizing one candidate-glucokinase (GCK). We show here that overexpression of GCK inhibits the proliferation of HCC cells via induction of intracellular lactate accumulation and subsequently causes energy crisis due to NAD+ depletion. This suggests GCK functions as a tumor suppressor gene and may be involved in HCC development. In conclusion, these data provide valuable clues for further investigations of the process of tumorigenesis in human cancer.

Reprogramming enriches for somatic cell clones with small-scale mutations in cancer-associated genes.

Cellular therapies based on induced pluripotent stem cells (iPSCs) come out of age and an increasing number of clinical trials applying iPSC-based transplants are ongoing or in preparation. Recent studies, however, demonstrated a high number of small-scale mutations in iPSCs. Although the mutational load in iPSCs seems to be largely derived from their parental cells, it is still unknown whether reprogramming may enrich for individual mutations that could lead to loss of functionality and tumor formation from iPSC derivatives. 30 hiPSC lines were analyzed by whole exome sequencing. High accuracy amplicon sequencing showed that all analyzed small-scale variants pre-existed in their parental cells and that individual mutations present in small subpopulations of parental cells become enriched among hiPSC clones during reprogramming. Among those, putatively actionable driver mutations affect genes related to cell-cycle control, cell death, and pluripotency and may confer a selective advantage during reprogramming. Finally, a short hairpin RNA (shRNA)-based experimental approach was applied to provide additional evidence for the individual impact of such genes on the reprogramming efficiency. In conclusion, we show that enriched mutations in curated onco- and tumor suppressor genes may account for an increased tumor risk and impact the clinical value of patient-derived hiPSCs.

Sequence diversity of the Pseudomonas aeruginosa population in loci that undergo microevolution in cystic fibrosis airways.

Five hundred and thirty-four unrelated Pseudomonas aeruginosa isolates from inanimate habitats, patients with cystic fibrosis (CF) and other human infections were sequenced in 19 genes that had been identified previously as the hot spots of genomic within-host evolution in serial isolates from 12 CF lungs. Amplicon sequencing confirmed a significantly higher sequence diversity of the 19 loci in P. aeruginosa isolates from CF patients compared to those from other habitats, but this overrepresentation was mainly due to the larger share of synonymous substitutions. Correspondingly, non-synonymous substitutions were either rare (gltT, lepA, ptsP) or benign (nuoL, fleR, pelF) in some loci. Other loci, however, showed an accumulation of non-neutral coding variants. Strains from the CF habitat were often mutated at evolutionarily conserved positions in the elements of stringent response (RelA, SpoT), LPS (PagL), polyamine transport (SpuE, SpuF) and alginate biosynthesis (AlgG, AlgU). The strongest skew towards the CF lung habitat was seen for amino acid sequence variants in AlgG that clustered in the carbohydrate-binding/sugar hydrolysis domain. The master regulators of quorum sensing lasR and rhlR were frequent targets for coding variants in isolates from chronic and acute human infections. Unique variants in lasR showed strong evidence of positive selection indicated by d N/d S values of ~4. The pelA gene that encodes a multidomain enzyme involved in both the formation and dispersion of Pel biofilms carried the highest number of single-nucleotide variants among the 19 genes and was the only gene with a higher frequency of missense mutations in P. aeruginosa strains from non-CF habitats than in isolates from CF airways. PelA protein variants are widely distributed in the P. aeruginosa population. In conclusion, coding variants in a subset of the examined loci are indeed characteristic for the adaptation of P. aeruginosa to the CF airways, but for other loci the elevated mutation rate is more indicative of infections in human habitats (lasR, rhlR) or global diversifying selection (pelA).

Novel candidate genes for ECT response prediction-a pilot study analyzing the DNA methylome of depressed patients receiving electroconvulsive therapy.

BACKGROUND: Major depressive disorder (MDD) represents a serious global health concern. The urge for efficient MDD treatment strategies is presently hindered by the incomplete knowledge of its underlying pathomechanism. Despite recent progress (highlighting both genetics and the environment, and thus DNA methylation, to be relevant for its development), 30-50% of MDD patients still fail to reach remission with standard treatment approaches. Electroconvulsive therapy (ECT) is one of the most powerful options for the treatment of pharmacoresistant depression; nevertheless, ECT remission rates barely reach 50% in large-scale naturalistic population-based studies. To optimize MDD treatment strategies and enable personalized medicine in the long- term, prospective indicators of ECT response are thus in great need. Because recent target-driven analyses revealed DNA methylation baseline differences between ECT responder groups, we analyzed the DNA methylome of depressed ECT patients using next-generation sequencing. In this pilot study, we did not only aim to find novel targets for ECT response prediction but also to get a deeper insight into its possible mechanism of action. RESULTS: Longitudinal DNA methylation analysis of peripheral blood mononuclear cells isolated from a cohort of treatment-resistant MDD patients (n = 12; time points: before and after 1st and last ECT, respectively) using a TruSeq-Methyl Capture EPIC Kit for library preparation, led to the following results: (1) The global DNA methylation differed neither between the four measured time points nor between ECT responders (n = 8) and non-responders (n = 4). (2) Analyzing the DNA methylation variance for every probe (=1476812 single CpG sites) revealed eight novel candidate genes to be implicated in ECT response (protein-coding genes: RNF175, RNF213, TBC1D14, TMC5, WSCD1; genes encoding for putative long non-coding RNA transcripts: AC018685.2, AC098617.1, CLCN3P1). (3) In addition, DNA methylation of two CpG sites (located within AQP10 and TRERF1) was found to change during the treatment course. CONCLUSIONS: We suggest ten novel candidate genes to be implicated in either ECT response or its possible mechanism. Because of the small sample size of our pilot study, our findings must be regarded as preliminary.

Intraclonal competitive fitness of longitudinal cystic fibrosis Pseudomonas aeruginosa airway isolates in liquid cultures.

The metabolically versatile Pseudomonas aeruginosa inhabits biotic and abiotic environments including the niche of cystic fibrosis (CF) airways. This study investigated how the adaptation to CF lungs affects the within-clone fitness of P. aeruginosa to grow and persist in liquid cultures in the presence of the clonal ancestors. Longitudinal clonal P. aeruginosa isolates that had been collected from 12 CF donors since the onset of colonization for up to 30 years was subjected to within-clone competition experiments. The relative quantities of individual strains were determined by marker-free amplicon sequencing of multiplex PCR products of strain-specific nucleotide sequence variants, a novel method that is generally applicable to studies in evolutionary genetics and microbial ecology with real-world strain collections. For 10 of the 12 examined patient courses, P. aeruginosa isolates of the first years of colonization grew faster in the presence of their clonal progeny than alone. Single growth of individual strains showed no temporal trend with colonization time, but in co-culture, the early isolates out-competed their clonal progeny. Irrespective of the genetic make-up of the clone and its genomic microevolution in CF lungs, the early isolates expressed fitness traits to win the within-clone competition that were absent in their progeny.

Effects of Lumacaftor-Ivacaftor Therapy on Cystic Fibrosis Transmembrane Conductance Regulator Function in Phe508del Homozygous Patients with Cystic Fibrosis.

RATIONALE: The combination of the CFTR (cystic fibrosis transmembrane conductance regulator) corrector lumacaftor with the potentiator ivacaftor has been approved for the treatment of patients with cystic fibrosis homozygous for the Phe508del CFTR mutation. The phase 3 trials examined clinical outcomes but did not evaluate CFTR function in patients. OBJECTIVES: To examine the effect of lumacaftor-ivacaftor on biomarkers of CFTR function in Phe508del homozygous patients with cystic fibrosis aged 12 years and older. METHODS: This prospective observational study assessed clinical outcomes including FEV1% predicted and body mass index, and CFTR biomarkers including sweat chloride concentration, nasal potential difference, and intestinal current measurement before and 8-16 weeks after initiation of lumacaftor-ivacaftor. MEASUREMENTS AND MAIN RESULTS: A total of 53 patients were enrolled in the study, and 52 patients had baseline and follow-up measurements. After initiation of lumacaftor-ivacaftor sweat chloride concentrations were reduced by 17.8 mmol/L (interquartile range [IQR], -25.9 to -6.1; P < 0.001), nasal potential difference showed partial rescue of CFTR function in nasal epithelia to a level of 10.2% (IQR, 0.0-26.1; P < 0.011), and intestinal current measurement showed functional improvement in rectal epithelia to a level of 17.7% of normal (IQR, 10.8-29.0; P < 0.001). All patients improved in at least one CFTR biomarker, but no correlations were found between CFTR biomarker responses and clinical outcomes. CONCLUSIONS: Lumacaftor-ivacaftor results in partial rescue of Phe508del CFTR function to levels comparable to the lower range of CFTR activity found in patients with residual function mutations. Functional improvement was detected even in the absence of short-term improvement of FEV1% predicted and body mass index. Clinical trial registered with www.clinicaltrials.gov (NCT02807415).

Impact of sample processing on human airways microbial metagenomes.

Whole metagenome shotgun sequencing provides information about the gene content and the composition of microbial communities provided that the processing of the samples does not introduce a methodology-driven bias. We tested the impact of DNA isolation and storage period on the metagenome profile. Deep throat swabs were collected from healthy adults and an infected infant. DNA was isolated by sonification or enzymatic lysis either immediately or after 24h storage in agar gel Amies transport medium at room temperature. Disruption of cells and subsequent fragmentation of DNA by sonification was as suitable as the common enzymatic lysis to generate high-quality metagenomes particularly for low total DNA input of less than ten nanograms. Conversely, storage of samples for 24h produced severely distorted metagenomes. The majority of species became less abundant or even extinct, whereas a few Streptococcus, Neisseria and Haemophilus spp. proliferated so that the total number of bacterial reads increased at the expense of human reads. We recommend that samples for metagenome analysis should be immediately processed or frozen at -80°C.

Multilocus amplicon sequencing of Pseudomonas aeruginosa cystic fibrosis airways isolates collected prior to and after early antipseudomonal chemotherapy.

BACKGROUND: Early antimicrobial chemotherapy can prevent or at least delay chronic cystic fibrosis (CF) airways infections with Pseudomonas aeruginosa. METHODS: During a 10-year study period P. aeruginosa was detected for the first time in 54 CF patients regularly seen at the CF centre Hannover. Amplicon sequencing of 34 loci of the P. aeruginosa core genome was performed in baseline and post-treatment isolates of the 15 CF patients who had remained P. aeruginosa - positive after the first round of antipseudomonal chemotherapy. RESULTS: Deep sequencing uncovered coexisting alternative nucleotides at in total 33 of 55,284 examined genome positions including six non-synonymous polymorphisms in the lasR gene, a key regulator of quorum sensing. After early treatment 42 of 50 novel nucleotide substitutions had emerged in exopolysaccharide biosynthesis, efflux pump and porin genes. CONCLUSIONS: Early treatment selects pathoadaptive mutations in P. aeruginosa that are typical for chronic infections of CF lungs.

The cystic fibrosis lower airways microbial metagenome.

Chronic airway infections determine most morbidity in people with cystic fibrosis (CF). Herein, we present unbiased quantitative data about the frequency and abundance of DNA viruses, archaea, bacteria, moulds and fungi in CF lower airways. Induced sputa were collected on several occasions from children, adolescents and adults with CF. Deep sputum metagenome sequencing identified, on average, approximately 10 DNA viruses or fungi and several hundred bacterial taxa. The metagenome of a CF patient was typically found to be made up of an individual signature of multiple, lowly abundant species superimposed by few disease-associated pathogens, such as Pseudomonas aeruginosa and Staphylococcus aureus, as major components. The host-associated signatures ranged from inconspicuous polymicrobial communities in healthy subjects to low-complexity microbiomes dominated by the typical CF pathogens in patients with advanced lung disease. The DNA virus community in CF lungs mainly consisted of phages and occasionally of human pathogens, such as adeno- and herpesviruses. The S. aureus and P. aeruginosa populations were composed of one major and numerous minor clone types. The rare clones constitute a low copy genetic resource that could rapidly expand as a response to habitat alterations, such as antimicrobial chemotherapy or invasion of novel microbes.

Molecular Epidemiology of Mutations in Antimicrobial Resistance Loci of Pseudomonas aeruginosa Isolates from Airways of Cystic Fibrosis Patients.

The chronic airway infections with Pseudomonas aeruginosa in people with cystic fibrosis (CF) are treated with aerosolized antibiotics, oral fluoroquinolones, and/or intravenous combination therapy with aminoglycosides and ß-lactam antibiotics. An international strain collection of 361 P. aeruginosa isolates from 258 CF patients seen at 30 CF clinics was examined for mutations in 17 antimicrobial susceptibility and resistance loci that had been identified as hot spots of mutation by genome sequencing of serial isolates from a single CF clinic. Combinatorial amplicon sequencing of pooled PCR products identified 1,112 sequence variants that were not present in the genomes of representative strains of the 20 most common clones of the global P. aeruginosa population. A high frequency of singular coding variants was seen in spuE, mexA, gyrA, rpoB, fusA1, mexZ, mexY, oprD, ampD, parR, parS, and envZ (amgS), reflecting the pressure upon P. aeruginosa in lungs of CF patients to generate novel protein variants. The proportion of nonneutral amino acid exchanges was high. Of the 17 loci, mexA, mexZ, and pagL were most frequently affected by independent stop mutations. Private and de novo mutations seem to play a pivotal role in the response of P. aeruginosa populations to the antimicrobial load and the individual CF host.

Intraclonal genome diversity of the major Pseudomonas aeruginosa clones C and PA14.

Bacterial populations differentiate at the subspecies level into clonal complexes. Intraclonal genome diversity was studied in 100 isolates of the two dominant Pseudomonas aeruginosa clones C and PA14 collected from the inanimate environment, acute and chronic infections. The core genome was highly conserved among clone members with a median pairwise within-clone single nucleotide sequence diversity of 8 × 10(-6) for clone C and 2 × 10(-5) for clone PA14. The composition of the accessory genome was, on the other hand, as variable within the clone as between unrelated clones. Each strain carried a large cargo of unique genes. The two dominant worldwide distributed P. aeruginosa clones combine an almost invariant core with the flexible gain and loss of genetic elements that spread by horizontal transfer.

Extended spectrum of human glucose-6-phosphatase catalytic subunit 3 deficiency: novel genotypes and phenotypic variability in severe congenital neutropenia.

OBJECTIVE: To delineate the phenotypic and molecular spectrum of patients with a syndromic variant of severe congenital neutropenia (SCN) due to mutations in the gene encoding glucose-6-phosphatase catalytic subunit 3 (G6PC3). STUDY DESIGN: Patients with syndromic SCN were characterized for associated malformations and referred to us for G6PC3 mutational analysis. RESULTS: In a cohort of 31 patients with syndromic SCN, we identified 16 patients with G6PC3 deficiency including 11 patients with novel biallelic mutations. We show that nonhematologic features of G6PC3 deficiency are good predictive indicators for mutations in G6PC3. Additionally, we demonstrate genetic variability in this disease and define novel features such as growth hormone deficiency, genital malformations, disrupted bone remodeling, and abnormalities of the integument. G6PC3 mutations may be associated with hydronephrosis or facial dysmorphism. The risk of transition to myelodysplastic syndrome/acute myeloid leukemia may be lower than in other genetically defined SCN subgroups. CONCLUSIONS: The phenotypic and molecular spectrum in G6PC3 deficiency is wider than previously appreciated. The risk of transition to myelodysplastic syndrome or acute myeloid leukemia may be lower in G6PC3 deficiency compared with other subgroups of SCN.